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Abstract 

O. 

^i^ ' A multidimensional study of local multiplicity fluctuations and multiparticle correlations 

Q^ ■ of hadrons produced in T? decays is performed. The study is based on the data sample 

^> \ of more than 4x10® events recorded with the OPAL detector at LEP. The fluctuations 

j><I \ and correlations are analysed in terms of the normalized scaled factorial moments and 

cumulants up to the fifth order. The moments are observed to have intermittency-like be- 
On. haviour, which is found to be more pronounced with increasing dimension. The large data 
(— I \ sample allows for the first time a study of the factorial cumulants in e"'"e~ annihilation. 
The analysis of the cumulants shows the existence of genuine multiparticle correlations 
with a strong intermittency rise up to higher orders. These correlations are found to be 
rS \ stronger in higher dimensions. The decomposition of the factorial moments into lower- 
ed ' order correlations shows that the dynamical fluctuations have important contributions 
from genuine many-particle correlations. The Monte Carlo models Jetset 7.4 and Her- 
WIG 5.9 are found to reproduce the trend of the measured moments and cumulants but 
they underestimate the magnitudes. The results are found to be consistent with QCD 
jet formation dynamics, although additional contributions from other mechanisms in the 
hadronization process cannot be excluded. 
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1 Introduction 

Particle density fluctuations of hadronic final states produced in high-energy collisions have 
been extensively investigated in the last decade. For recent reviews, see e.g. Refs. |1|, 0. 
Dynamical (i.e., non-statistical) fluctuations were observed, establishing the phenomenon of 
intermittency, the increase of factorial moments with decreasing bin size p|. The intermit- 
tency approach of studying the distributions of particles in restricted regions of phase space 
allows a detailed analysis of the dynamics of hadropro duct ion. Furthermore, the behaviour 
of the factorial moments shows the self-similar nature of density fluctuations, i.e., the parti- 
cle distributions show similar fluctuations on all resolution scales, a characteristic of fractals 0]. 

Despite numerous experimental and theoretical studies, the origin of intermittency remains 
unclear, although important features of this effect have been observed 0. For example, ex- 
perimental investigations have shown an enhancement of the phenomenon in e"'"e~ annihilation 
as compared to hadronic and nuclear collisions. Furthermore, larger intermittency effects have 
been observed when several dynamical variables are considered together, as compared to the 
effect seen in one-dimensional analyses. Existing Monte Carlo (MC) models, which use parton 
shower simulations and differing fragmentation and hadronization models, simulate most de- 
tails of e"'"e~ collisions and the general properties of hadronic interactions well, but fall short in 
predicting the intermittent structure found in the data, both in e"'"e~ collisions and in hadronic 
interactions. Theoretical approaches have not clarified the origin of the observed dynamical 
fluctuations. The intermittent behaviour of particle distributions may prove to be a strong 
test of QCD, which already provides guidelines for explaining the "soft" character of intermit- 
tency |,§. 

The goal of this study is to investigate the dynamical correlations of many-particle systems 
produced in e"'"e~ annihilation. One must be careful to separate out the effects of lower-order 
correlations when searching for higher-order ones. For example, a correlation in the production 
of pairs of particles in neighboring regions of phase space will necessarily induce correlations 
when particles are considered three at a time. It has been suggested that intermittency should 
therefore be analysed in terms of the factorial cumulant moments to reveal "genuine" multi- 
particle correlations by not being sensitive to the contributions of lower-order correlations 0. 
The investigations carried out for heavy-ion reactions have not shown any correlations higher 
than two-particle ones, while in studies of hadron-hadron collisions significant higher-order cor- 
relations have been observed, although the latter have been seen to weaken with increasing 
multiplicity at a fixed centre-of-mass energy 0. These effects have been explained as a conse- 
quence of the events consisting of superpositions of multiple independent particle sources p|. 
These findings suggest that interactions with a low number of very hard scattering processes, 
such as high-energy e"'"e~ annihilation, might be more sensitive to genuine multiparticle corre- 
lations. 

This paper describes the study of the intermittency phenomena and the genuine multiparti- 
cle correlations of charged particles in the three-dimensional phase space of rapidity, transverse 
momentum, and azimuthal angle, as defined in Section ^. This analysis uses more than four mil- 
lion multihadronic events recorded by the OPAL detector at the LEP collider with ^/s ~ mzo. 



This data sample is much larger than that used in OPAL's previous publication |^, and in 
other e"'"e~ investigations carried out at the Z° peak |1I0|-[14| and at lower energies [[I3l , p!5| . The 



statistical precision of this data sample allows us to extend the former intermittency analysis |^ 
to a multidimensional one with the possibility to reach high-order fluctuations at very small 
bins. With this high statistics we are able for the first time to search for genuine multiparticle 
correlations by means of normalized factorial cumulants in e"'"e~ collisions. 

The paper is organised as follows. In Section the normalized factorial moments and 
factorial cumulant moments are introduced and the technique of extracting dynamical fluctua- 
tions and genuine multiparticle correlations is given. The detector, data sample and correction 
procedure are described in Section |. The results and their comparison with MC predictions 
are presented and discussed in Section |[ In Section ^ both a summary and conclusions are 
presented. 

2 The method 

2.1 Scaled factorial moments and intermittency 

To search for local dynamical fluctuations we use the normalized scaled factorial moments 
introduced in Ref. [Q. The moments are defined as 

F, = Ar'^(;5)/iVM. (1) 

Here Um is the gth order factorial multinomial, nm{nm — 1) ■ ■ • (^m — Q' + 1), with rim par- 
ticles in the mth bin of the phase space (e.g. rapidity interval) divided into M equal bins. 
Nm is the number of particles in the mth bin summed over all the M events in the sample, 
^m = J2j=ii^m)j- The bar indicates averaging over the bins in each event, (1/M) X]m=i 
("horizontal" averaging), while the angle brackets denote averaging over the events ("vertical" 
averaging) . 

The moments in Eq. (|l]) are given in the modified form, in contrast to those used in the 
earlier e"'"e~ studies [P-[T5[|. This form has been suggested in Ref. [T^ to take into account the 



bias arising from the assumption of infinite statistics in the normalization calculation P, |T7 
and affects the moments, particularly those computed with small bins. 



These moments, defined in Eq. (|I|), are the so-called "horizontal" moments |17|, where the 



fluctuations are first averaged over all M bins in each event and then the average over all events 
is takenj^ These moments are determined for a flat-shape single-particle distribution. In order 
to account for non- flatness, we apply the correction procedure proposed in Refs. []T^,|T^, so that 



the corrected modified factorial moments 16] are given in the reduced form. 



F^ = FJR,, R, = Nt^/Nl:i\ (2) 



^ For a survey of types of factorial moments see Ref. [p^ 



where Rg is the correction factor, and it is equal to unity for a flat single-particle distribution. 

The non-statistical fluctuations have been shown to lead to increasing factorial moments 
with increasing number (decreasing size) of the phase-space regions, or bins. Such an increase, 
expressed as a scaling law, 

Fg{M) oc M^« (M -^ oo), 0<(pg^q-l, (3) 

indicates the presence of self-similar dynamics. This increase is called intermittency and the 
powers (pq, or intermittency slopes, show the strength of the effect. The size of the smallest 
bin under investigation is limited by the characteristic correlation length (saturation effect) 
and/or by the apparatus resolution ||^. In practice, saturation happens much earlier because of 
statistical limitations (the "empty bin effect" [^), which we here attempt to reduce by using 



the modified moments. Therefore, in the following we use the term "intermittency" to refer 
only to the rise of the factorial moments with decreasing bin size. 



2.2 Factorial cumulant moments and genuine multiparticle correla- 
tions 

To extract the genuine multiparticle correlations, the technique of normalized factorial cumu- 
lant moments, or cumulants. Kg, proposed in Ref. |^, is used. The cumulants are constructed 
from the g-particle cumulant correlation functions which vanish whenever one of their argu- 
ments becomes independent of the others, so that they measure the genuine correlations. The 
factorial cumulants remove the influence of the statistical component of the correlations in the 
same way as the factorial moments. 

Similarly to the factorial moments, we use the corrected modified cumulants, defined as 



K^=Af'^k^"^^/Nt\ (4) 



The normalization factor, A',^ (instead of Nm), comes from the correction procedure expressed 
in Eq. (§) and takes into account the non-flat shape of the single-particle distribution. The 



kg factors are the unnormalized factorial cumulant moments or the Mueller moments 121 



and represent genuine g-particle correlations where the lower-order contributions are eliminated 
by subtracting the appropriate combinations of the unnormalized factorial moments, (rim), of 
order p < q from the gth order one, e.g. 

kt^ = {n^3-3{n^^){n^)+2{n^f. (5) 

In order to find contributions from genuine multiparticle correlations to the factorial mo- 



ments we use the relations between the moments and the cumulants ||22|| , 



F2 = K2 + 1, 

Fs = K3 + 3K2 + I, (6) 

F4 = 7^4 + 47^3 + 3(K^'"^)2 + 6i^2 + 1, 



'M Tv-M I 1 n ?<-. I 1 c; ,' i^{™) 



i^5 + 5/^4 + 10 K^'Kr' + IOK3 + 15 {KrT + WK2 + 1, etc. 



'M _ i,!"*) 



with i^g = fcq /{nm}'^- Here and below in this section, for brevity we omit the superscript 
C. These relations are exact for the "vertical" moments and cumulants and the errors intro- 
duced by using them with corrected horizontal moments and cumulants are negligibly small, 
except when very high orders are considered for variables whose single-particle distributions 
are markedly non- uniform 0. 

The composition of the g-particle dynamical fluctuations from the genuine lower-order p- 
particle correlations is tested by the comparison of the gth order factorial moments with the 
p-particle contribution Fq calculated by the above Eqs. @, which are truncated up to the 
i^p-terms. The excess of Fg over Fq demonstrates the importance of correlations of order 
higher than p in the measured g-particle fluctuations. 

For example, the two-particle contribution F^ to the fourth-order factorial moment can be 
expressed as 

fP =3{Kt^^ + 6K2 + l, (7) 

while the contribution from two and three-particle correlations, F^ , is 



-(2+3) _ Ai^ , Q/7^M 



Fr > = 47^3 + 3(irrO' + QK2 + 1. (8) 



3 Experimental details 



The data used in this study were recorded with the OPAL detector |2^ at the LEP e"'"e~ collider 
at CERN. The analysis is restricted to charged particles measured in the central tracking cham- 
bers. The inner vertex detector has a high precision in impact parameter reconstruction. The 
large diameter jet chamber and outer layer of longitudinal drift chambers allow accurate mea- 
surements in the planes perpendicular and parallel to the beam axis. The jet chamber provides 
up to 159 space points per track, and allows particle identification by measuring the ionisa- 
tion energy loss, dE/dx, of charged particles |23|. All tracking chambers are surrounded by a 



solenoidal coil providing a magnetic field of 0.435 T along the beam axis. The resolution of the 
component of momentum perpendicular to the beam axis is (j{pt)/pt = A/(0.0023pi)^ + (0.018)^ 
for pt in GeV/c, and the resolution in the angle 9 between the charged particle's direction and 
the electron beam is a{9) = b mrad within the acceptance of the analysis presented here. In 
multihadronic events, the ionisation energy loss measurement has been obtained with a resolu- 
tion of ~ 3.5% for tracks with 159 measured points. 



The present study was performed with a sample of approximately 4.1x10^ hadronic Z° de- 
cays collected from 1991 through 1995. About 96% of this sample was collected at the Z° peak 
energy and the remaining part was collected within ±3 GeV of the peak. Over 98% of charged 
hadrons were detected. 

The event selection criteria are based on the multihadronic event selection algorithms de- 
scribed in Refs. [0,|2^, and are similar to those used in other LEP intermittency studies [pD|-p^. 



For each event, "good" charged tracks were accepted if they 

• had at least 20 measured points in the jet chamber; 

• had a first measured point closer than 70 cm from the beam axis; 

• pass within 5 cm of the e"'"e~ collision point in the projection perpendicular to the beam 
axis, with the corresponding distance along the beam axis not exceeding 40 cm; 

• had a momentum component transverse to the beam direction greater than 0.15 GeV/c; 

• had a momentum smaller than 10 GeV/c; 

• had a measured polar angle of the track with respect to the beam direction satisfying 
Icos^l < 0.93; 

• had a mean energy loss, dE/dx, in the jet chamber smaller than 9 keV/cm to reject 
electrons and positrons. 

Selected multihadron events were required to have 

• at least 5 good tracks; 

• a momentum imbalance of < 0.4:^/s, which is defined as the magnitude of the vector sum 
of the momenta of all charged particles; 

• the sum of the energies of all tracks (assumed to be pions) greater than 0.2^/s^, 

• I cos6's| < 0.7, where 9s is the polar angle of the event sphericity axis with respect to the 
beam direction. The sphericity axis is calculated using all good tracks and electromagnetic 
and hadronic calorimeter clusters. 

The first three requirements provide rejection of background from non-hadronic Z° decays, 
two-photon events, beam-wall interactions, and beam-gas scattering. The last requirement en- 
sures that the event is well contained in the most sensitive volume of the detector. A total of 
about 2.3x10^ events were selected for further analysis. 

We also used two samples of about 2x10^ simulated events each, generated at Z° peak 
energy using the following MC generators: 



Jetset version 7.4 [26| with the parton shower followed by string formation and frag- 
mentation, 



Herwig version 5.9 |2^ with the parton shower followed by cluster fragmentation. 



The parameters of both MC codes were tuned with OPAL data p8[ to provide a good descrip- 



tion of the distributions of the measured event-shape variables and single-particle kinematic 
variables. 

In this study we chose rapidity, azimuthal angle and transverse momentum as individual 
track kinematic variables. These are frequently used in multihadronic studies |2^ and, in par- 
ticular, for intermittency and correlation analyses |jl|,^^,|,0. To make our study compatible 



with other investigations carried out in e"'"e annihilations, these variables are calculated with 
respect to the sphericity axis as follows: 

• Rapidity, y = 0.5\n[{E + p\^)/{E — P||)] with E and p\\ being the energy (assuming the 
pion mass) and longitudinal momentum of the particle in the interval —2.0^y^2.0. 

• Transverse momentum in the interval — 2.4^1n(pT)^0.7. The log scale is introduced, since 
the exponential shape of the pT-distribution causes instability in the average multiplicity 
calculations, even for pt bins of intermediate size. 

• Azimuthal angle, $, calculated with respect to the eigenvector of the momentum tensor 
having the smallest eigenvalue, in the plane perpendicular to the sphericity axis. The 
angle $ is defined in the interval 0^$^27r. 

The single-particle distributions of the data sample and of the MC (corrected to the hadron 
level, see below) are shown in Fig. |l|. In the following study the maximum number of bins is 
taken to be Mmax=400, so that the one-dimensional minimal bin size of the above kinematic 
variables are: 5|/min=0.01, 5$min=0.9° and (5(lnpT)mm— 0.008. M^ax is the same as was chosen 
in our previous intermittency publication [§ and is the largest value of M used so far. The 
experimental resolution of each variable was estimated by a MC simulation. It was found that 
the OPAL detector allows the study of an intermittency signal down to distances of the above 
mentioned minimal bin sizes, although detector effects become important for bin sizes less than 
0.04 in rapidity, smaller than 3° in azimuthal angle and less than 0.02 in Inp^. The distribu- 
tions in several dimensions have the advantage that the event phase space may be subdivided 
into many more bins than M^ax before the limits of detector resolution are reached. As in our 
former analysis, the smallest bin sizes used are found not to affect the observations. 

To correct the measured moments for the effects of geometrical acceptance, kinematic cuts, 
initial-state radiation, resolution and particle decays, we apply the correction procedure adopted 
in our earlier factorial moment study 0| and in analogous investigations done with e^e~ annihi- 
lation [p!0|-p!^ . Two samples of multihadronic events were generated using the Jetset 7.4 MC 
program. The first sample does not include the effects of initial-state radiation, and all parti- 
cles with lifetimes longer than 3xl0~^° seconds were regarded as stable. The generator-level 
factorial moments and cumulants are calculated directly from the charged particle distributions 
of this sample without any selection criteria. The second sample was generated including the 

9 



effects of finite lifetimes and initial-state radiation and was passed through a full simulation of 



the OPAL detector |31I]. The corresponding detector- level moments were calculated from this 
set using the same reconstruction and selection algorithms as used for the measured data. The 
corrected moments were then determined by multiplying the measured ones by the factor 

U,{M) = X^{M)g,jX,{M),,t , (9) 

with Xg = F^ or K^ defined by Eqs. (|) and (§. 

The correction factors Uq have also been computed using the Herwig event sample. For 
both Jetset and Herwig MC generators, the correction factor tends to be less than unity 
and rises with order q as has been observed earlier ^,^,[1^. The difference between the Xg 



quantities calculated with the Jetset and Herwig generated samples has been used in the 
estimation of the systematic uncertainties. The statistical uncertainties on the Jetset Uq fac- 
tor have been also incorporated into the systematic uncertainties in this analysis. 

Another contribution to the systematic uncertainties has been evaluated by changing the 
above track and event selection criteria. The moments have been calculated from the data 
sample of about two million events with the following variations in the selection criteria: the 
first measured point was required to be closer than 40 cm to the beam axis, the requirement 
of the transverse momentum with respect to the beam axis was removed, the total momentum 
was required to be less than 40 GeV/c, the charged track polar angle acceptance was changed 
to I cos^l < 0.7, and the requirement on the mean energy loss was removed. The deviations of 
the moments with these changes modify the results by no more than a few percent and do not 
influence their behaviour. 

The total errors have been calculated by adding the systematic and statistical uncertainties 
in quadrature. The systematic uncertainties are shown separately in the figures (except in those 



given in Section [4.3| ) and are dominant at large M. Statistical uncertainties based on the MC 



samples are similar to those obtained from the data. 

It was verified that the results do not appreciably change if one removes from the analysis 
those events which were taken at energies off the Z^ peak energy. 

4 Results and discussion 

4.1 The measurements 

In this section we present the factorial moments F^ and the cumulants Kg , defined in Eqs. (0) 
and (^, respectively, and compute them in the yx^xpT phase space and its projections. The 
moments are shown in Figs. to ^ and the cumulants are given in Figs. |^ to as a function of 
M, the number of bins in one, two and three dimensions. Both the factorial moments and the 
cumulants are measured up to the fifth order. The second-order cumulants are not shown since 
their behaviour is determined directly by the second-order factorial moments, as can be seen 
from Eqs. (^). The higher-order cumulants behave differently from the same-order moments 
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because the cumulants contain combinations of lower-order fluctuations, which are taken into 
account in Eq. (^ by means of the Mueller moments as in Eq. (^). 

Overall, the factorial moments and the cumulants depend very similarly on the bin width. 
The factorial moments obey the scaling-law of Eq. (^) in almost all the phase-space projections 
(except in the p^-subspace), and the cumulants show analogous intermittent behaviour up to 
high orders. This behaviour becomes more pronounced when the analysis is extended to several 
dimensions where, in contrast to the one- dimensional case, no saturation with decreasing bin 
size is observed. The largest moments and cumulants and the largest intermittency slopes are 
found in the ?/x$ subspace. 

The saturation at small bin sizes observed in the one-dimensional analysis in rapidity and 
azimuthal angle (Figs. || and |) agrees with that predicted by QCD ^ and is a consequence 
of the transition to the regime where the running of the QCD coupling as comes into play. 
The dynamics governing particle density fluctuations in small bins occurs at low energy scales 
with larger values of as- One can see that the moments and the cumulants actually have steep 
linear rises at M<20 (5|/>0.2, (5<l>>18°) and level off at large values of M. Figs. ^ and ^ show 
that the transition point shifts to larger M (smaller bin sizes) as the order q of the moments 
increases, also in accordance with QCD calculations. 

The enhancement of the intermittency effect in higher dimensions as shown in Figs. ^ and 
^ for the factorial moments and in Figs. |] and ^ for the cumulants, has been attributed to 



"shadowing" i.e., studies in lower dimensions lose information due to projection [^, ^. A 
model of emission of strongly coUimated particles, clustered in both rapidity and azimuth, has 
been suggested in [^. In the framework of this so-called "pencil-jet" model, a strong increase 
of the factorial moments is expected in the |/x$ subspace compared to y or $ separately. Al- 
though formation of such jet structures has been confirmed experimentally, the increase has 
been found to be much less than that predicted [Q]. 

Jet structure also explains the behaviour observed in the y^px and (^xpx plots of Figs. |] 
and|^. Indeed, the moments and the cumulants in yxpr and ^xp^ for the same total M are not 
found to increase as compared to those in y and $, respectively, since there is no intermittency 
in the transverse momentum subspace. Similarly, the moments in ?/x$ are approximately equal 
to those in yx^xp^ at the same total M. At the same time, the intermittency is seen for larger 
M in higher dimensions, indicating the presence of the dynamical fluctuations and correlations 
in additional dimensions. 



In Ref. I^y] it was claimed that the increase of the factorial moments with the addition of 
new dimensions is a trivial consequence of a phase-space factor and has nothing to do with the 
jet formation mechanism. Our observations show that this is true if one compares the moments 
at the same abscissa M^^'^, where d is the dimension of the subspace. However, comparing 
multidimensional moments (and the cumulants) at the same total number of bins, one obtains 
contributions of self-similar structure from different projections, as is the case with the jet- 
structure contributions to the ?/x$ moments. 
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The values of the cumulants are positive in most of the cases, indicating that multiparticle 
dynamical correlations indeed are significantly present in the particle-production process. Large 
values of the cumulants of the order g^4 are seen in ?/x$ and yx^xprp. Non-zero high-order 
cumulants are also found in rapidity and yxpT- This shows that the factorial moments at ^=5, 
the highest order considered here, have important contributions from lower-order correlations. 



a point discussed in Section 4.3 



4.2 Comparison with the Monte Carlo models 

In Figs. H through |^ the data are compared with the predictions of the Jetset and Herwig 
MC models. Both MC models describe the general behaviour of the factorial moments and 
cumulants and show significant positive multiparticle correlations. 

The one-dimensional factorial moments (Fig. H) and cumulants (Fig. §) in rapidity and in 
azimuthal angle show that while the MC describe the data rather well at small M (large bin 
sizes), the models tend to fall below the data, starting at intermediate M. The discrepancy 
rises with M and with the order of the moments q. The saturation effect sets in earlier in the 
MC models than it does in the data. In the transverse momentum projection, the models show 
quite different behaviour. Jetset underestimates the moments and the cumulants, whereas 
Herwig strongly overestimates them. Figs. |^, ^ ^ and |^ show that there is a better agreement 
between the data and the MC in high dimensions. 

From these comparisons one can conclude that both MC models reproduce the data well 
while neither of them is particularly preferred. The perturbative parton shower, on which 
both MC models are based, seems to play an important role in the origin of the dynamical 
fluctuations and correlations in e"'"e~ annihilation. The observed differences between the two 
MC descriptions indicate that the last steps of the hadronization process are not described 
correctly 0. Contributions from additional mechanisms to the observed fluctuations and cor- 
relations are not excluded. 

4.3 Contributions from multiparticle correlations 

This section describes the contributions of genuine multiparticle correlations to dynamical fluc- 
tuations. To this end we compare in Figs. |8|-[12| the measured corrected factorial moments to 



the lower-order contributions, Fq , calculated using Eqs. (|^). 

Fig. 1^ shows the one-dimensional case. The fluctuations in transverse momentum are not 
shown since they do not exhibit intermittency behaviour (see Fig. |^). A significant contribu- 
tion of high-order genuine correlations appears. Two-particle correlations in rapidity and in 
azimuthal angle are insufficient to explain the measured three-particle fluctuations. At g=4, 
four-particle correlations are also necessary. The importance of four-particle genuine correla- 
tions is also demonstrated by the flve-particle fluctuations in the rapidity subspace, where the 
addition of the fourth-order cumulants becomes essential. The flfth-order moment study cannot 
be performed for the azimuthal angle variable because the non-uniformity of the $ spectrum 
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leads to large values of the correction factor Rr, which makes the relations (6) inapplicable. The 
difference between the moments and the correlation contributions increases with decreasing bin 
size. 



The genuine multiparticle contributions, also seen to be important in the two-dimensional 
?/x$ analysis, are shown in Fig. ^ The failure of the genuine two- and three-particle corre- 
lations to describe the four-particle dynamical fluctuations indicates a significant four-particle 
contribution in the observed high-order fluctuations. The need to account for higher-order 
correlations is visible also in F^ for small bin-sizes. The comparison of F^ and F^ for other 



two-dimensional subspaces, yxpr and ^xpx, is shown in Figs. |Ty and [11] and also indicates a 
considerable contribution from three-particle genuine correlations. As in the one- dimensional 
case in $, the large i?„ factor for g>3 does not allow use of Eqs. (BI) in these cases. The same 



is seen in the three-dimensional study (Fig. |T2|) for g=5, whereas for q < 5, the contribution of 
multiparticle correlations up to the fourth order is well illustrated. 



4.4 Comparison with other experiments 

4.4.1 Factorial moments in e+e annihilation 

Studies of intermittency in e+e^ interactions have been carried out mainly in one-dimension 
using the rapidity variable P,p!0|-[T5[]. The first three-dimensional analysis of factorial moments 
was perfomed by CELLO [0 in Lorentz-invariant phase space. DELPHI |l^ has presented a 
three-dimensional analysis of intermittency in yx^xpT phase space and its projections, as it 
is chosen in this paper. The values of the moments and their M- dependence, found in all these 
studies, are similar to those obtained here. 



In one-dimensional rapidity and azimuthal angle the saturation of the factorial moments 
has been observed at the same M in all investigations |T^-|I^. The moments in two and three 



dimensions have been found [T^ to be larger and to have steeper intermittency slopes than 
those in one dimension. Jet evolution has been suggested [^,T^, 3^| as the main source of mul- 
tiparticle fluctuations, similarly to our finding. However, the saturation shown by DELPHI for 
the two-dimensional factorial moments are not present in our analysis due to the high statistics 
and the modification used to take into account the contents of small bins. 



In agreement with recent e"'"e studies |10-12, 14|, the MC models used here are found to 



describe the behaviour of the measured moments, although they underestimate their magni- 
tudes. 

4.4.2 Factorial moments and cumulants in hadronic collisons 

The factorial moments have also been investigated in yx^xpT phase space and in its projec- 



tions in hadron-hadron collisons by NA22 BB 



In the NA22 pj-subspace analysis, hadronic interactions have shown a visible intermittency 
effect, in contrast to its absence in e^e~ annihilation. In $, on the other hand, one finds sensi- 
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tive dynamical fluctuations in e"'"e~ collisions, while in hadronic collisions the fluctuations are 
strongly suppressed by a statistical component. No saturation has been observed in y and pt 
subspaces in hadronic collisions. 

In several dimensions the factorial moments computed in hadronic interactions also differ 
from those in e"'"e~ annihialtion. In two-dimensions, the largest moments are found to be in 
yxpT, and the intermittency effect is observed to be the strongest in yx^. These moments 
show a faster increase (larger ipg) than those in one dimension, although their values are lower 
and saturate already at q=3. The three-dimensional moments show strong increase with de- 
creasing bin size, although their values are found to be closer to the two-dimensional moments. 
In contrast to the linear log-log plots of the scaling behaviour (§) observed in e"'"e~ annihilation, 
three-dimensional factorial moments in hadronic collisions scale only approximately; they rise 
more quickly than the power law of Eq. (j^). This difference is attributed to the difference 
in dynamics between soft and hard processes that leads to isotropic dynamical fluctuations in 
e"'"e~ annihilation and to anisotropic ones in hadron-hadron collisions ||36||. 



Large cumulant values are expected in e"'"e~ annihilation due to a small number of very 
energetic particle sources, as was mentioned in the introduction. Indeed, the cumulants mea- 
sured here are much higher than those found in hadronic interactions. Furthermore, they have 
non-zero values for higher orders, while in hadronic collisions they are consistent with zero at 
g>3. Also, the power-law increase in M is seen to be stronger in the present data. 



In hadronic interactions 0, both in one [Q and several dimensions |^, the factorial 
moments were found to be basically composed of two-particle correlations. Our study shows 
sensitive contributions of correlations of orders g=3 and even 4 to the dynamical fluctuations in 
e"'"e~ annihilation. The observed contributions increase with decreasing bin size, in agreement 
with trends seen in hadronic interactions. 

5 Summary and Conclusions 

A multidimensional study of local fluctuations and genuine multiparticle correlations in hadronic 
decays of the Z° is carried out with the four-million event sample of the OPAL data collected 
at the LEP collider. The sphericity axis is chosen as a reference axis, and the phase space 
is deflned by the rapidity, azimuthal angle and transverse momentum variables. The analysis 
is based on the intermittency approach and represents an investigation of the normalized fac- 
torial moments and, for the flrst time in e"'"e~ annihilation, the normalized cumulants. The 
quantities studied have been corrected to reduce the bias due to the non-uniform shapes of the 
single-particle distributions and have been modifled to eliminate effects of finite statistics. The 
factorial moments and cumulants are computed up to the fifth order and down to very small 
bin sizes. 

The factorial moments show an intermittency behaviour which is strongly enhanced as the 
dimension of the subspace increases from one to three. In the one-dimensional analysis, the 
intermittency signal is found to be larger in rapidity than in azimuthal angle and to vanish in 
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transverse momentum. The moments in rapidity and in azimuthal angle saturate at interme- 
diate bin sizes, in agreement with the QCD expectation for the transition to the running Ug 
couphng regime. No saturation is observed in two and three dimensions, a consequence of jet 
formation. 

Our study of the cumulants shows that they have large positive values, indicating the ex- 
istence of genuine multiparticle correlations. The cumulants are found to be much larger than 
those in hadronic interactions, suggesting an increase of the correlations with the decrease of 
the number of independent subprocesses present in the reaction. The cumulants are analysed 
in subspaces of different dimensions, and the genuine multiparticle correlations are found to 
be larger in higher dimensions. This study reveals genuine correlations up to four particles 
in one dimension and significant five-particle correlations in higher dimensions. Large higher- 
order correlations measured in the |/x<l> two-dimensional projection confirm the jet structure 
of dense groups of particles. The cumulants show intermittency rises which are stronger than 
for the corresponding factorial moments. 

The large statistics of the present analysis allows the observation of contributions of many- 
particle correlations to the measured dynamical fluctuations. Considerable contributions up to 
four-particle are observed in the expansion of the factorial moments in terms of cumulants. The 
contributions increase with decreasing bin size reflecting the underlying self-similar dynamics. 

The measurements are compared to Herwig 5.9 and Jetset 7.4 predictions. In general, 
these Monte Carlo models are found to reproduce the behaviour of the moments and the cumu- 
lants, while underestimating the measured values starting at intermediate bin sizes. High-order 
multiparticle correlations are found to be present in both models used. The observations confirm 
jet structure formation as an important contribution to the correlations, but other mechanisms 
of the hadronization process are possibly also relevant. 
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Figure 1: The single-particle rapidity, the logarithm of the transverse momentum, and azimuthal 
angle distributions in the data and as predicted in two Monte Carlo models. All kinematic 
variables are calculated with respect to the sphericity axis. The distributions are corrected for 
detector effects in a bin-by-bin manner. The arrows show the intervals used in the analysis. 
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Figure 2: Factorial moments of order q = 2 to 5 vs. the number of bins M of the one- 
dimensional rapidity, azimuthal angle, and transverse momentum subspaces in comparison with 
the predictions of two Monte Carlo models. The total error is shown along with the statistical 
error (inner error bars) for each point. 
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Figure 3: Factorial moments of order q = 2 to 5 vs. vM, where M is the number of bins 
of the two-dimensional subspaces of rapidity, azimuthal angle, and transverse momentum com- 
binations, in comparison with the predictions of two Monte Carlo models. The total error is 
shown along with the statistical error (inner error bars) for each point. 



21 



OPAL 



"i^ 10 ' 



10 



10 



10 



yxOxp. 



Data 

JETSET 7.4 
HER WIG 5.9 



} 




10 



10 



M 



1/3 



Figure 4: Factorial moments of order q = 2 to 5 vs. \fM , where M is the number of bins 
of the three-dimensional rapidity, azimuthal angle, and transverse momentum phase space, in 
comparison with the predictions of two Monte Carlo models. The total error is shown along 
with the statistical error (inner error bars) for each point. 
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Figure 5: Cumulants of order q = 3 to 5 vs. the number of bins M of the one- dimensional 
rapidity, azimuthal angle, and transverse momentum subspaces in comparison with the predic- 
tions of two Monte Carlo models. The total error is shown along with the statistical error (inner 
error bars) for each point. 
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Figure 6: Cumulants of order q = 3 to 5 vs. y M , where M is the number of bins of the two- 
dimensional subspaces of rapidity, azimuthal angle, and transverse momentum combinations, 
in comparison with the predictions of two Monte Carlo models. The total error is shown along 
with the statistical error (inner error bars) for each point. 
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Figure 7: Cumulants of order q = 3 to 5 vs. -vM where M is the number of bins of the three- 
dimensional rapidity, azimuthal angle, and transverse momentum phase space, in comparison 
with the predictions of two Monte Carlo models. The total error is shown along with the 
statistical error (inner error bars) for each point. 
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Figure 8: Decomposition of the factorial moments F^ into multiparticle correlation contribu- 
tions Fq in the one- dimensional rapidity and azimuthal angle subspaces. The error bars show 
the total uncertainties. 
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Figure 9: Decomposition of the factorial moments F^ into multiparticle correlation contribu- 
tions Fq in the two-dimensional subspace of rapidity and azimuthal angle. The error bars 
show the total uncertainties. 
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Figure 10: Decomposition of the factorial moments F^" into multiparticle correlation contribu- 
tions Fq in the two-dimensional subspace of rapidity and transverse momentum. The error 
bars show the total uncertainties. 
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Figure 11: Decomposition of the factorial moments F^ into multiparticle correlation contribu- 
tions Fq in the two-dimensional subspace of azimuthal angle and transverse momentum. The 
error bars show the total uncertainties. 
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Figure 12: Decomposition of the factorial moments F^ into multiparticle correlation contribu- 
tions Fq in the three dimensions. The error bars show the total uncertainties. 
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